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Abstract The limitations for the coherent manipulation of 
neutral atoms with fabricated solid state devices, so-called 
'atom chips', are addressed. Specifically, we examine the 
dominant decoherence mechanism, which is due to the mag- 
netic noise originating from the surface of the atom chip. It is 
shown that the contribution of fluctuations in the chip wires 
at the shot noise level is not negligible. We estimate the co- 
herence times and discuss ways to increase them. Our main 
conclusion is that future advances should allow for coherence 
times as long as 1 second, a few /im away from the surface. 

PACS: 03.75. -b Matter waves - 32.80.Lg Mechanical effects 
of light on atoms and ions - 03.67.Lx Quantum computation 
- 05. 40. -a Fluctuation phenomena and noise 



In the quest for physical implementations of quantum infor- 
mation processing, "atom chips" are currently of great inter- 
est. This is because they promise well-controlled quantum 
optical manipulations of neutral atoms in integrated and scal- 
able microtrap arrays. In these traps, atoms are strongly con- 
fined by electromagnetic fields close to nanostructured solid- 
state substrates. Microtraps used in current experiments are 
magnetic traps produced by current-carrying wires [jj l |2| , 3|||, 
|[|[0,|J^,[Io[ and periodically magnetized substrates ipl |, or 
hybrid traps involving optical or electric fields Jl2|, |l3l| . 

In this paper, we discuss the limitations that wire-based 
magnetic traps on atom chips may have to face when they are 
downscaled into the micron range. Recently, both theoretical 
and experimental indications have been found that the 'hot' 
chip substrate — typically held at room temperature — is not 
a quiet environment: at distances below a few 100 /im from 
the chip, the trap lifetime is shorter than in free space and the 
atom temperature increases [||^||l4[|l5,16|. While it is not 
excluded that strong compression in these microtraps plays 
a role due to enhanced collisional interactions (see 0171 ] for 
a review), noisy magnetic fields may also be involved. They 
provide a coupling to the environment that may cause loss, 
heating and decoherence and are elaborated upon in this pa- 
per. We review the sources of magnetic fields and quote esti- 
mates for trapping and coherence times. In particular, we dis- 



cuss the contribution of electric current noise at the shot noise 
level and evaluate its spectral and spatial properties. This is 
compared to the noise due to the thermal chip substrate. 



1 Atom chip 'building block': the side guide 

In the 1930's, Frisch and Segre realized that when a homoge- 
neous magnetic field ('bias field') is superimposed with the 
field of a straight wire current, the magnetic field vanishes on 
a line parallel to the current (see figure [j]) [|l8|]. In the vicinity 
of this line, the field increases in a quadrupolar fashion. The 
height of the field zero is given by 



2tt B b ' 



(1) 



where = 47rmmG/A is the vacuum permeability, I the 
wire current and Bb the bias field magnitude. (This expres- 
sion applies to finite size wires provided their diameter <C h.) 
This 'side guide' can be implemented on an atom chip using 
a lithographically etched wire on the chip surface. The setup 
is ideal for miniaturization since the guide height is reduced 
using smaller currents (with less ohmic dissipation). 

Atoms with a nonzero magnetic moment fi are trapped 
around the magnetic field minimum provided fi keeps an 
antiparallel orientation with respect to the local magnetic 
field (adiabatic approximation). These 'weak field seekers' 
are attracted to the field zero and get trapped in a potential, 
/i|||Bo(x)|, proportional to the magnitude of the magnetic 
field. (The subscript distinguishes the static trapping field 
from the field fluctuations discussed later.) It is particularly 
interesting to introduce a longitudinal magnetic field (along 
the wire) so that the field magnitude never reaches zero: this 
helps satisfying the adiabatic approximation everywhere and 
reduces so-called 'Majorana flips' into untrapped magnetic 
sublevels. The trapping potential is then a harmonic well in- 
stead of a linear one around the guide center. This is the guide 
geometry we focus on in the following. The relevant param- 
eters are the trap height h (Eq.[l]), the oscillation frequency 
fl/2-K in the harmonic well, and the Larmor frequency cjl = 
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jtii||Bo(r)|/fi where B (r) is the field at the trap center and 
fx\\ is the magnetic moment along the center field [p[|lO|]. For 
the lowest quantum states in the potential well, the magnetic 
field is predominantly longitudinal, and the single-particle 
wave functions are approximately harmonic oscillator states 
for quantum numbers up to n ~ uj^/ Q. 

The validity of the adiabatic approximation is also deter- 
mined by the ratio Lo^j Q: it has been estimated that nonadi- 
abatic losses are exponentially suppressed when this ratio is 
large (see also []2(J). Similarly, tunnelling losses down 
to the chip surface can be made exponentially small with a 
sufficiently high and 'thick' potential barrier. The main loss 
channel is then provided by the uncontrollable coupling to the 
surface via magnetic noise. 



2 Interaction with magnetic noise 

2.1 Noise spectrum and spin flip rate 

The coupling of the atomic magnetic moment to fluctuating 
magnetic fields gives rise to both spin flips and changes in 
the center-of-mass motion (scattering). The rate of these pro- 
cesses is given by Fermi's Golden Rule. We recall here that it 
can be conveniently expressed in terms of the noise spectrum 
of the magnetic field fluctuations. (See [|T]| for a similar ap- 
proach and Chap. IV of [ p2| ] for the derivation of a full master 
equation.) 

If we write |i) and |f) for the atomic states before and 
after the transition, the transition rate is 

Jl-f = y 5>(I) I (F, f \H int |I, i) | 2 S(E F + Et-Ek- Ei), 

F,I 

(2) 

where |I) and |F) are initial and final states for the field, the 
summation being an average over the initial field states (with 
probabilities p(l)) and a trace over the final field states. The 
interaction Hamiltonian is given by H m t = fi ■ B(x). 

Consider first the rate for spin flips. Since only a subset 
of magnetic sublevels |mi) are weak field seekers, spin flips 
— > \rri[ ) are responsible for trap loss. The magnetic field 
is evaluated at the position r of the trap center. (An average 
over the atomic position distribution would be more accu- 
rate.) We write the (5-function for energy conservation as a 
time integral over e KEi-E F -huj n )t/h w jjere htun = E f - Ei. 
The exponential e 1 ( Bl ~ BF ) i / s can be removed by introducing 
the field operators in the Heisenberg picture 

B(r, t) = e i ^*/ ?i B(r)e- iHot / ?i (3) 

(here, Hq is the free field Hamiltonian) and taking matrix ele- 
ments of this operator between the initial and final field states. 
This gives 

j&-z*)t/*Q\B(r)\F) = (I|B(r,t)|F). (4) 



The sum over the final states |F) now reduces to a complete- 
ness relation and we get (a, (3 denote field components) 

2TrhJ2p( l )( l \B a (r)\F)(F\B fj (r)\l)S(E F - E 1 ~ huj) 

F,I 

oo 

= f &tj" t Y J pm\Bo t {r,t)B {T,m 

= S a p(r;io). (5) 

In the last line, we have defined the magnetic noise spectrum 
which is the Fourier transform of the field's autocorrelation 
function. The rate for spin flips can now be written as 

Zi_>f = -j (m i |/i Q |mf)(mf|/z / 3|mi)S , Q ^(r, -w fi ). 

OL.j3 — x.y ,Z 

(6) 

Since nif ^ mu the matrix elements of p a are only nonzero 
for directions perpendicular to the magnetic field at the trap 
center. We also recover the selection rule nif — rrii — ± 1 so 
that the relevant transition frequency is the Larmor frequency 
l^fil = ^l- The spin flip rate gives the order of magnitude of 
trap loss even if more than one weak-field seeking Zeeman 
states, m\ = +2, +1, say, are trapped (possible with many 
of the alkali atoms). This is because the matrix elements be- 
tween adjacent sublevels do not significantly differ in magni- 
tude so that the atoms reach the non-trapped sublevel rri{ = 
after a time ~ 2/1^2^+1- 

We finally note that as long as the behaviour of the 'envi- 
ronment' (the field) is ignored in the description of the atom's 
dynamics, the noise spectrum is the only quantity needed 
to characterize the environment. It is also an experimentally 
measurable quantity: for example, the rms magnetic noise 
(-B^(r)) 1 / 2 measured by a spectrum analyzer in a given fre- 
quency band Alo/2tt around uj is (2S xx (r, uj) Aj/271") 1 / 2 , 
the factor 2 accounting for the sum over positive and nega- 
tive frequencies. The atomic spin flip rate may be regarded as 
an alternative way to measure the noise spectrum. In order of 
magnitude, the magnetic moment is comparable to the Bohr 
magneton, p-Q (pb/^ttTI — 1.4MHz/G), and we get 

r^ t (r)~oms-\»/» B f Sa ^\ (7) 

pi /Hz 

Note that current SQUID magnetometers are able to detect 
magnetic field noise even on the 10 iT/\/Hz scale [p3|]. 

2.2 Near field noise 

The magnetic noise spectrum close to a solid substrate shows 
dramatic differences with respect to the well-known black- 
body spectrum. This is because thermally excited currents in 
the substrate (Johnson-Nyquist noise) generate electromag- 
netic fields with substantial nonpropagating components. As 
an illustration, figure || shows the noise spectrum at a fixed 
distance (h = 1 p,m) above a metallic half-space, calculated 
along the lines of [O]. One gets an increase of several orders 
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of magnitude for low frequencies (wavelength A much larger 
than h). In addition, the spectrum is rather flat in this range. 
Only at high frequencies (A < h), the Planck spectrum is re- 
covered because the detector then enters the far field of the 
source. 

In the low-frequency range, the magnetic near field noise 
above a planar substrate is approximately characterized by 
the spectrum [0] 



— «1: 
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(8) 

where T, g are the substrate temperature and resistivity, re- 
spectively, h is the observation distance, s a p = diag(^, 1, |) 
is a diagonal tensor (the distinguished axis is along the sur- 
face normal), and 8{u) = y / 2g//j,oUj is the skin depth. In 
order of magnitude, the corresponding spin flip rate is sur- 
prisingly large for traps at a micrometer distance: 



h < S(lo l ) 



100 s 



_ 1 W MB ) 2 (T/300K) 
(g/gcu)(h/nm) 



(9) 



where gc u — 1-7 x 10~ 6 f2cm is the copper resistivity. 
The flip rate is plotted as a function of distance in fig- 
ure ^|. A larger Larmor frequency (longitudinal bias field) 
only helps reducing spin flips when h > <5(wl)> giving a 
scaling 2i_+f ocw, ' . Finally, there are many plausible rea- 
sons to expect that the linear dependence on T of the flip 
rate (^) does not continue down to very low temperatures: 
the full Bose-Einstein occupation number has to be used, 
k^T i ► Tilji, (1 — e -fiuL/k B T^-i ^ an( j omer sources of mag- 
netic fields may come into play like spin waves, lattice vi- 
brations etc. One thus probably gets lifetimes (and coherence 
times, see the following) shorter than predicted by (Q). 

At room temperature, the estimate (Q) describes a worst 
case because most materials have a resistivity larger than cop- 
per. A structured substrate like a thin metallic layer or a wire 
also helps: the magnetic noise then decreases more rapidly 
with distance h [ |To| , [24| ]. As a general rule, the smaller the 
amount of metallic material, the lower the magnetic noise. 
This can be understood from a simple model where ther- 
mally excited currents in each volume element of the con- 
ducting substrate contribute to the total magnetic noise above 
(a similar approach has been used in [ |23| , |2~5| ]). For distances 
h <C 8(u>), one obtains a noise spectrum accurate within a 
factor of two from an incoherent addition of magnetostatic 
fields, neglecting the influence of the material on the field 
propagation. For example, at a distance h from a single, 
thin wire with radius a <C h, the theory of Ref. [Q] yields 
a magnetic noise tensor given by Eq.(||) with the replace- 
ment s a p/h i — ► (7ra 2 //i 3 )diag(2, i, |), showing a similar 
weak anisotropy. The components specify the azimuthal, ra- 
dial and longitudinal directions in cylindrical coordinates, re- 
spectively. The correction involving the skin depth 5(lo) can- 
not be obtained in the magnetostatic approximation. 

Experimental data [|[^,|l6|] indicate that trap lifetimes as 
long as or even longer than estimated by (^) (several 100 ms 
to 100 s) are achievable: these traps were built close to thin 



wires (3 to 90 /im width) or semiconductor substrates cov- 
ered with thin metallic layers (thickness <C h), at distances 
between 20 /im and 2 mm. It should be noted that the dis- 
tance dependence of the trap lifetime measured in JI^ ] only 
agrees qualitatively with our theory, and more detailed inves- 
tigations are required. 



2.3 Current noise 

The electric currents that generate the side guide fields are 
also subject to fluctuations that drive spin flips and deform the 
trapping fields. The impact of technical noise can be reduced 
using electronic filtering, "quiet" drivers, by correlating the 
currents for the wire and the bias field etc. This works down to 
the most fundamental level where the noise due to the discrete 
value of the electron charge comes into play ('shot noise'). 
For a current / in free space, shot noise has a frequency- 
independent spectrum given by 



SNi = el « 0.16 nA 2 /Hz 



(10) 



where e is the charge quantum. Note that currents in a solid 
wire can have fluctuations below the shot noise limit because 
the Coulomb interaction correlates the electrons. In the fol- 
lowing, we use SNj as a convenient reference value. 

If the wire current in a side guide has a noise spectrum 
Si(u), it creates a magnetic field with 



/x|el Sj(cj) 
4ir 2 h 2 SN! ' 



(ID 



This scalar quantity gives one component of the noise tensor. 
The spin flip rate is then of the order of 



Is 



{h/iim) 2 SNi A' 



(12) 



We conclude that miniaturization beyond the micrometer 
scale requires extremely low current noise to achieve trap 
lifetimes longer than seconds. It may well turn out that 
very large scale atom chip integration is only possible with 
static magnetic fields, generated by magnetized nanostruc- 
tures pT||. Experimental trap lifetime data [^.Q show that 
currently used power supplies are quiet enough not to reduce 
lifetimes at distances above ~ 10 /im. 



3 Spin coherence 

Coherent manipulation on atom chips requires both that 
atoms stay trapped and that their quantum state be preserved. 
Here, we focus on the influence of magnetic noise on the in- 
ternal (spin) states. For example, different magnetic sublevels 
or hyperfine states are interesting candidates to implement a 
bit of quantum information (qubit). Noise-induced transitions 
between sublevels erase the qubit, and this occurs on the same 
timescale as the spin flips discussed in the previous section. 
But the information contained in quantum superpositions can 
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also be lost by pure 'dephasing', without changing sublevel 
populations [26|. For this process, which opens an additional 
channel for decoherence, longitudinal magnetic fields, i.e. po- 
larized along the static trapping field ('phase noise'), come 
into play. 

Consider two magnetic sublevels |mi) and \1ri2) that are 
simultaneously trapped. A magnetic field fluctuation parallel 
to the static trapping field changes the energy difference be- 
tween the two qubit states by 



AEQt) = An\\B\\(T,t) 



(13) 



where B\\ is the longitudinal magnetic field and A/j,n = 
(m2\n\\ \m2) — ( TO i|M||l m i) tne differential magnetic mo- 
ment. (This difference can be substantially reduced in alka- 
lis by choosing hyperfine states that only differ in the nuclear 
spin state.) Let us now expose a superposition of \m\) and 
1 7712) to pure magnetic phase noise during an interaction time 
t. The off-diagonal element of the corresponding 2x2 density 
matrix is then proportional to [^6[] 



exp — / At' AE(t') ) = exp 



Afih 



(14) 



4 Vibrational coherence 

In this section, we turn to scattering processes that leave the 
atoms in the magnetic trap, but perturb their center-of-mass 
motion. This occurs whenever the magnetic noise is not spa- 
tially homogeneous. A simple model is suggested showing 
that the typical length scale for magnetic inhomogeneities is 
of the order of the trap height h for both near field noise and 
current noise at the shot noise level. We discuss the relation 
between random changes in the atoms' momentum and the 
decoherence of their density matrix in position space. 



4.1 Scattering rate 

We again use Fermi's Golden Rule (||), but now the initial 
and final states are given by wave functions ^i f(x). In the 
following, we retain only a single trapped magnetic sublevel 
|mf) = |toi) and assume that the magnetic moment has the 
same orientation for all relevant center-of-mass states. Its ma- 
trix elements then reduce to m (the component along the 
trapping field). Writing out the overlap integral between wave 
functions, the transition rate (0) becomes 



where S\\ (r; 0) is the low-frequency limit of the magnetic 
noise spectrum. (More precisely, one needs the spectrum av- 
eraged over the frequency range . . . 1/t. We neglect this 
complication since the relevant spectra are flat in this range.) 

From Eq.(jT4j), we conclude that dephasing leads to expo- 
nential decoherence of qubit superpositions with a rate simi- 
lar to the spin flip rate (cf. Eq.(f5()). The decoherence rate gets 
smaller when the logical states have the same magnetic mo- 
ment (reducing ziyu.11), or when the magnetic field component 
£>H shows much less noise. 

Near field noise is rather isotropic, as shown by the spec- 
trum (|j), and therefore contributes equally to spin flips and 
phase noise. Noise in the wire current gives only fluctua- 
tions perpendicular to the guide axis so that dephasing is sup- 
pressed close to the guide center. Note that this suppression 
is not complete because of the finite, transverse width of the 
trapped wave function. In addition, gravity can displace the 
actual trap center with respect to the magnetic field mini- 
mum (the 'gravitational sag' familiar from Bose condensates 
in magnetic traps, see also [^[]). In the harmonic approxima- 
tion, a simple calculation leads to a reduction of the dephas- 
ing rate by a factor [M g//j,\\b) 2 where in order of magnitude 



Mg 01 (M/amu)[g/(10m/s 2 )] 
/i||6 ' (/*||//i B )[6/(G/cm)] 



(15) 



This ratio can be made quite small (< 10~ 6 ) using typical 
magnetic gradients b = B^/h achievable with the side guide. 
Magnetic near fields thus remain as the main source of de- 
phasing noise, with a rate basically scaling like the spin flip 
rate. 



2 

Li r 

r^i = /d 3 xdVA^(x)M fi (x')S||(x,x';- Wfi ), 



where the wavefunction overlap is given by 



M fi (x)=#(x)Vi(x), 



(16) 



(17) 



the transition energy hiva is the difference between the cen- 
ter of mass levels, and the field correlation spectrum is the 
generalization of the noise spectrum (Bh: 



00 

5||(x,x»= f die i - t ^p(I)(I|S||(x,i)i?||(x',0)|I). 
J 1 



(18) 

A useful figure that can be extracted from this function is 
the correlation length l c that governs the variation of ( |l8| ) as 
a function of distance s = x — x'. As a general rule, one 
has Sii(x, x!;oj) — > if s 3> l c because the fields Bn(x,t) 
and £?||(x',t) become decorrelated. More quantitatively, we 
define here l c by the following expansion for small deviations 
from the trap center r 



S\\ (x, x'; u) « S\\ (r;w) 



(x- 



(19) 



provided |x — r|, |x' — r| <C l c 



In Sec.^, we show that the correlation length l c is comparable 
to the trap height h, for both near field and current noise. 
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4.2 Heating 

Consider the vibrational motion along one transverse direc- 
tion (the x-axis, say) in the harmonic region of the side 
guide. Magnetic noise can induce transitions between dif- 
ferent quantum states in the trap ('heating'). The transition 
— * 1 between the ground and first excited states is particu- 
larly interesting and involves the overlap integral (u>n = ft) 

2 

dxdx' Mf i (x)M n (x')S l \{x,x';-n) « ^-S||(r; -ft) 

(20) 

where a — (K/2Mi2) 1 ' 2 is the (rms) size of the trap ground 
state (M is the atomic mass) and the noise spectrum is eval- 
uated at the trap center r. Eq.(pO|) is derived in the limit 
of strong confinement (a much smaller than the correlation 
length l c ), using the expansion (|l9|). 

Comparing ( p^ ) to (^), we conclude that heating in tight 
traps is suppressed relative to spin flips, due to the small ratio 
(a/l c ) 2 <C 1. We stress that this is not due to the increase in 
the trap frequency (because the noise is essentially white), but 
due to the small size of the trap ground state. Since l c ~ h, the 
excitation rate Jb-^i follows power laws 1 /h 3 , 1 /h 4 , instead 
of the flip rates (||, |l2|), depending on the noise source. For 
example, the rate due to near field noise above a metallic half- 
space scales like 



r 



Is" 1 (/V/z B ) 2 (T s /300K) 



(M/amu)(J7/27rl00kHz)( £ »/ £ -cu)(V^ m ) 3 



(21) 

Heating is also relevant for the decoherence of qubits imple- 
mented in different vibrational levels. One can derive a mas- 
ter equation for the density matrix in the harmonic well that 
shows that the corresponding off-diagonal elements decay at 
a rate comparable to Jb-^i. For details, see jl5[|2^ ]. 

Finally, we note that even spatially homogeneous mag- 
netic fluctuations can induce heating when they change the 
trap position or curvature. Displacing the trap is equivalent to 
a force and drives the transition — > 1 with a rate [ |To| , |2"7| , 



r ^ = ^s h (-n) 



(22) 



where M is the atomic mass and Sh(oj) is the spectrum of 
the trap height fluctuations. Taking into account only field 
fluctuations caused by the wire current, we get in order of 
magnitude 



r ^i - 3s- 1 (Af/amu)(ft/27rl00kHz) 3 
I/A Si{Q) 
X (B b /G) 2 SNj • 



(23) 



This rate is still reasonably small for sufficiently large bias 
fields B\, > 50 G. Note, however, that a very strong confine- 
ment may not be possible due to the increase with the trap fre- 
quency ft: this occurs because a larger spring constant Mft 2 
translates position fluctuations into larger forces. 

Fluctuations in the trap position due to technical noise can 
be reduced by correlating the currents in the wire with those 



procuding the bias fields. The trap curvature, however, then 
still fluctuates and changes the oscillation frequency ft. This 
gives a parametric resonance on the 0^2 transition with a 
rate plplES]: 



r. 



i 



0^2 



Sn(-2f2) 



3xir8rl (^100kHz) 2 S;(2fi) 



I/A 



SN! 



(24) 



(25) 



which is substantially smaller than the rate (£3|). Fluctuations 
in ft also induce phase noise on the quantum states in the 
harmonic trap because their energy difference is TiQ. Arguing 
as in Sec.^], one finds a dephasing rate comparable to Zo->2. 

In experiments with trapped atoms, it is relatively sim- 
ple to observe the rate of temperature increase, T. In the 
harmonic approximation, one finds a 'heating rate' fc^T = 
ftftlo^i for noise driving the transitions n — ► n ± 1 p7[|2^]. 
Currently observed values are in the range of 0.05— 1 /iK/s [|5[ 
|l6| l . This is orders of magnitude larger than predicted by ( plf ) 
and excludes an origin dominated by near field noise. Heat- 
ing can probably be attributed to current fluctuations, as given 
by (|23"|), assuming typical power supply noise spectra (above 
shot noise). Ambient electromagnetic noise ('electrosmog') 
may also play a role and is currently under investigation. 



5 Spatial coherence 

5.7 Scattering 'cross section' 

We now consider the quasi-free motion along the side guide 
axis Oz, the transverse motion is assumed to be 'frozen out'. 
Noise induces scattering p; — > pf between different mo- 
mentum states, where the wavefunction overlap Mn(z) = 
L _1 exp(— iqnz) involves the wavevector transfer Kqn = 
Pf—Pi {L is a normalization length). The transition frequency 
is ujfi — qapi/M + hq 2 i /2M. We end up with the spatial 
Fourier transform of the correlation function (the vectors x, 
x' only differ in their z-components) 



V\\ fdzdz' . 



h 2 J L 2 

,2 



.e^-^Siifcx's-wfi) 



« ^^(rj-wfi) j^-e^ s C{s--u f {). (26) 

To perform the last step, we have assumed that the correla- 
tions involve only the distance s = z — z' (statistically homo- 
geneous noise) and used the normalized correlation function 

«7(0;w) = l): 



C(z-z';uj) = 



S||(x,x';a;) 
5||(r;w) • 



(27) 



Finally, we get a transition rate per wavevector transfer 
dr/dq by dividing by the wavevector spacing dq = 2n/L 
in the quantization volume: 



dn 



dq 



(28) 
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where the rate 7 is the prefactor of the integral in Eq.(pq). 
If the elements of the magnetic noise tensor are of compara- 
ble magnitude, this rate is of the same order as the spin flip 
rate (||). As mentioned at the end of Sec.^J, this is the case for 
near field noise and to a some extent also for fluctuations due 
to current noise. 

We conclude from the differential scattering rate ( |28| ) that 
typical momenta exchanged with the noise field have a mag- 
nitude qa ~ h/l c given by the inverse noise correlation 
length. 



5.2 Decoherence 

The relation to spatial decoherence of the atoms has been 
made more quantitative in [^40, where the following master 
equation for the density matrix W(z, p) in the Wigner repre- 
sentation is derived 



d, 



p_ 

M' 



d z )W(z,p) 



dq 



dq~ 



(W(z,p + hq) - W(z,p) 



(29) 



The scattering integral on the right hand side describes pro- 
cesses p p + fiq, while the left hand side gives the free 
ballistic motion along the guide axis. 

The master equation ( 29 ) can be solved exactly using the 
fact that the 'scattering cross section' ( ^ ) is essentially inde- 
pendent of the transition frequency u)a for the relevant mag- 
netic field fluctuations. One finds the following expression for 
the spatially averaged coherence function of the atoms [E4|| 



p(s,t) = J dz(i/j*(z + s,t)ip(z,t)) (30) 
= p(s,0)exp(- 7 *[l-C(s;0)]), (31) 

where the brackets denote the average with respect to the 
noise and C(s; 0) is the low-frequency limit of the noise cor- 
relation function. 

From the coherence function (|3l|), we identify j[l — 
C(s; 0)] as the decoherence rate for spatially separated su- 
perposition states: for a splitting greater than the correlation 
length, s ^S> l c , the correlation function C(s; 0) is zero and 
the superposition decays into a statistical mixture on a time 
scale given by 1/7, comparable to the spin lifetime. Super- 
positions with smaller splitting s decay more slowly, with a 
rate scaling like [1 — C(s;Q)] r y (s/Z c ) 2 7 -C 7, using the 
expansion (19). This behaviour was also found in a decoher- 
ence model by W. H. Zurek p9[ ] who used a master equation 
in Fokker-Planck form instead of our Eq.(p^). 

We note that the decoherence rate just found also de- 
scribes the dephasing between the arms of a guided matter 
wave interferometer Jl0|, ^o|pll ]. This follows from an argu- 
ment similar to that used in Eq.(|l4j). In this context, s is the 
(transverse) separation between the arms of the interferome- 
ter. For more details, we refer to JTofa . 



5.3 Decoherence of a condensate 

We now discuss the extension of the previous results to the 
case of a Bose condensate in a linear waveguide (see [^] 
for a review of experiments). The single-particle wave func- 
tions then have to be replaced by collective modes, and 
the effective potential is changed due to atom-atom inter- 
actions. For example, the (transverse) ground state has a 
larger width compared to the single-particle wave function so 
that the atomic spins are no longer aligned along the guide 
axis. We present here preliminary results for the decoher- 
ence of a Bose condensate focussing on a quasi- ID regime 
[H||H§||]. Adopting a mean-field description, we have 
performed Monte-Carlo simulations of the one-dimensional 
Gross-Pitaevski equation 



m 



dip 



h 2 ^ 



2M dz 2 



+ V(z,t)iP + g\iP(z)\ 2 iP (32) 



where the coupling constant g = 2hf2a s is proportional to the 
(transverse) trap frequency and the s-wave scattering length 
a s > (repulsive interactions). The random potential V(z, t) 
is chosen in accordance with the correlation functions rel- 
evant for atom chip traps (white noise and lorentzian spa- 
tial correlations). The initial situation is a condensate in the 
ground state of a harmonic trap superimposed on the wave- 
guide potential. Quantum (phase) fluctuations [[34j[35]] are ig- 
nored assuming effectively T at = 0. The harmonic confine- 
ment is instantaneously released at t = 0, and the cloud ex- 
pands along the waveguide axis. 

The simulation results given in figure [| show the spatially 
averaged coherence function of the condensate, Eq.(|30|), as a 
function of the separation s for different expansion times t. 
One observes that in the presence of noise (scattering rate 
7 7^ 0), the coherence length (the width of p(s, t)) is reduced 
as time increases — the cloud breaks apart in mutually inco- 
herent patches that have 'seen' different noise potentials. The 
dotted lines in the upper left figure give the prediction ( |3~l| ) of 
the master equation for noninteracting atoms, and we note a 
very good agreement between the analytics and the numerical 
data. 

We find that the decoherence scenario is not qualitatively 
changed by a moderate self-interaction, as a comparison of 
the panels for g — and g — 10 shows. Only for short times 
is the coherence length of an interacting cloud larger because 
the ground state is broadened by the interactions. A more de- 
tailed investigation of condensate decoherence, including an- 
alytical approximations and the limit of strong interactions, 
will be presented elsewhere [|(|. 



6 Noise correlation length 

We finally show that magnetic noise involved in atom chips 
has a correlation length l c ~ h, for both near field and current 
noise. 
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6.1 Near field noise 

The spatial correlation function for the electric near field 
above a thermal, planar substrate was studied in [|37|]. A simi- 
lar calculation yields for the magnetic near field the following 
normalized correlation function: 



C(a;0) 



8h 2 



16/i 2 /3 



(2h + Vs 2 +4h 2 )Vs 2 +'ih 2 s 2 + 16/i 2 /3 

(33) 

This gives the correlation for the magnetic field component 
Bn, taken at positions z, z' on the side guide axis with a sep- 
aration s. We have assumed that the field wavelength (fixed 
by the transition frequency ojh) is much larger than the rele- 
vant distances s, h. The lorentzian form in Eq.(p3|) is a good 
approximation for all distances where C(s;0) is sensibly 
nonzero and shows even more explicitly that the correlation 
length is of the order of the guide height h. This is because 
the noise fields radiated by each volume element of the sub- 
strate are quasi-static in the near field and decay algebraically 
(no retardation). 



6.2 Shot noise correlations 

One might think at first sight that magnetic noise due to cur- 
rent fluctuations should have a large correlation length be- 
cause the relevant electromagnetic frequencies (kHz to MHz 
range) propagate with a large wavelength along the wire. The 
analysis of near field fluctuations has shown, however, that 
the wavelength is not really the relevant scale as soon as 
one is sensitive to non-propagating fields generated by nearby 
sources. For this reason, we suggest a simple toy model for 
the flow of electrons through a thin wire that allows to recover 
both the noise spectrum of the magnetic field and its spatial 
correlations. 

The ingredients of the model are sketched in figure ^] and 
details are given in Appendix ^[ The electrons are assumed 
to move independently and ballistically (the Drude electron 
gas model), their transverse position in the wire is neglected 
compared to the guide distance h. It is beyond the scope 
of this model to describe correlations between the electrons 
that could lead to lower current fluctuations. Neither does the 
model describe diffusive electron transport, we comment on 
that below. 

The result of the model is the magnetic noise spec- 
trum ( |4l| ) given in the Appendix. We recover the previous 
spectrum (^lj) (with current noise at the shot noise level, 
Si(u) — SNj) in the low-frequency limit where u>h -C v 
for all relevant electron velocities v. At high frequencies, 
the noise is reduced because one needs fast electrons to pro- 
duce magnetic field 'pulses' with short duration ~ h/v. This 
is visible in Figure g where the spectrum (0), for fixed 
h = 1 im\, is shown as a function of frequency, normalized 
to its low-frequency limit. The high-frequency cutoff occurs 
at uih ^> wf where vp is the Fermi velocity when the Fermi- 
Dirac distribution is taken for P(v). A Maxwellian distribu- 
tion gives a similar behavior with a lower cutoff, as shown 



in the figure. Note again that at most frequencies relevant for 
atom chip traps, the noise spectrum can be assumed flat. 

The characteristic length scale l c for the spatial correla- 
tion of the shot noise fields can also be read off from Eq.(p4]). 
In the directions perpendicular to the guide axis, it is given 
by the guide height h, and we recover the same correlation 
length as for near field fluctuations. For the motion along the 
guide axis, the simplest way is compute the differential scat- 
tering rate (|28j). For a process p- } — > p; + hq, we get 

dP- 



-7SN-^Pfe/M) [qhK^qhtf 



(34) 



dq Mq 

where we have neglected the recoil shift fiq 2 /2M compared 
to the Doppler shift qp-jM. The scattering rate is 

Mjjcos 2 a ^2 eJ 



7SN 



(35) 



h 2 4tt 2 ^ 2 ' 

where a is the angle between the atomic spin and the wire 
field. If the harmonic approximation for the transverse mo- 
tion is not valid, cos a is not small and ( |35| ) is comparable 
to the spin flip rate (cf. Eq.|TTj), typically a few 1 s _1 . K\ in 
Eq.(p4|) is the modified Bessel function of the second kind, 
and P(v) is the electrons' velocity distribution, taken at the 
atomic velocity p-jM. It is interesting that the scattering in- 
volves a class of electrons co-moving with the atom: this sug- 
gests that the atomic wave in the linear guide is diffracted by 
the spatially confined field pulse. The divergence of the cross- 
section ( p4| ) for forward scattering (q — ► 0) is related to the 
long range behavior of the field pulse; this behavior also oc- 
curs for the Coulomb potential where the Rutherford cross 
section diverges in the forward direction. 

We can now conclude that also the longitudinal scattering 
is limited to momentum transfers hq < h/h. This is due to 
the large-argument asymptotics of the Bessel function in the 
result (B4l) 



qh»l: qhK^qh) ^ \j ^e~ qh , 



Tiqh 



(36) 



giving an exponential suppression for large qh. The shot 
noise field is thus also 'rough' on a scale l c ~ h, and does not 
behave qualitatively different compared to the thermal mag- 
netic near field. 

We finally comment on our neglect of the diffusive elec- 
tron motion in the wire. Since this effect gives rise to the 
nonzero resistivity of the wire, it is in fact included in the 
Johnson noise approach for the metallic substrate. The finite 
drift velocity of the electrons, which makes the distribution 
P(v) asymmetric, does not change our conclusions either be- 
cause it is typically much smaller than the width of P(v). The 
electron drift is taken into account in the spectrum shown in 
figure U 

7 Perspectives 

We have reviewed in this paper loss, heating and decoher- 
ence mechanisms for wire-based atom chips and their scal- 
ing with the microtrap geometry and the substrate material 
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properties. The importance of the shot noise level for current 
noise has been highlighted. Using a simple model, we have 
shown that the spatial correlation length of magnetic fields 
due to shot noise is fixed by the distance between microtrap 
and chip wire. 

The extreme miniaturization of atom chip traps below the 
1 scale may not be possible with "conventional", conduct- 
ing nanostructures, because magnetic field fluctuations due 
to thermal and technical current fluctuations become quite 
strong. The key process is a noise-induced change of the 
atomic sublevel \m) ('spin flip'). Its rate, which can be re- 
lated to the trap lifetime, also dictates the order of magnitude 
of more subtle processes involving heating of the center-of- 
mass motion or qubit dephasing. The timescale for useful co- 
herent manipulations is thus at least limited to a few 100 ms at 
a height of a few /im. Depending on the gate time of 2-qubit 
operations, this time scale may be sufficient. 

Several strategies leading to more robust atom chips can 
be imagined. Cooled substrates reduce thermal near field 
noise, with a gain in lifetime inversely proportional to the 
temperature. Quiet current drivers and/or superconducting 
wires are an option to reduce current noise, possibly below 
the shot noise level. Substrates with a permanent magnetiza- 
tion may also provide the required low-noise environment. 

The theory of magnetic noise close to complex, magne- 
tized structures can be developed starting from a description 
of the material in terms of its electric and magnetic suscepti- 
bilities. According to the fluctuation-dissipation theorem [23], 
[38| ], the imaginary parts of these fix the magnitude of the 
noise current and magnetization that generate thermal noise 
fields. This scheme can be used at low frequencies (high tem- 
peratures), which is the typical situation in atom chip traps, 
but can also be extended to high frequencies where the noise 
reduces to vacuum fluctuations, modified by the boundary 
conditions set by the substrate (see e.g. [ |39| ] for a review). 
This kind of approach can be used to compute atom-atom in- 
teractions mediated by virtual photon exchange. 

Another theoretical task is to estimate the coupling be- 
tween higher atomic levels and the electric and magnetic 
fields originating from the surface. For example, it has been 
proposed to make use of the electric dipole-dipole interaction 
to realize controlled two-qubit gates on atom chip traps p0| 
Rjj ]. The gate operation is sped up when the atoms are excited 
to high-lying Rydberg states, but these states also interact 
strongly with the chip substrate due to their large dipole mo- 
ments, hence the need for a review of their coherence times. 

Finally, the understanding of the impact of the surface on 
the fidelity of qubit operations needs to be studied in order to 
optimize the construction of an atom chip quantum processor. 
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A Correlation of magnetic fields generated by shot noise 

Each electron (charge e), during its passage with velocity 
v below the trapped atom, gives an electromagnetic 'pulse' 
whose vector potential is given by 



A(x, t) = e z e /c(x, t — to; v) 

_ e z /Joey/Air 



[x 2 + y 2 + (z - v(t - t )) 2 ] 1/2 ' 



(37) 



where the coordinates are chosen as shown in figure ||. The 
atom actually experiences an average vector potential that is 
due to the flow of many electrons, passing below the atom 
at random instants to- With the assumption of a stationary 
electron flow, we find 



<A(x,t)> 



y 2 ) + const. 



(38) 



where the average current is / = e (n) / At with n the num- 
ber of electrons flowing during the interval At. This average 
vector potential is time-independent and gives the static mag- 
netic field generated by the current. Note that it does not in- 
volve the velocity distribution of the electrons. (The calcula- 
tion is analogous to the calculation of a photodetector current, 
as outlined in chapter 9.8 of [0].) 

We are interested in the correlation function of the vector 
potential (from which the magnetic field correlations follow 
via differentiation). Subtracting the average value and per- 
forming again the average over the flowing electrons, we find 
using the approximation of independent electrons 

(A z (x,t + r)A z (x.',t)) - (A z (x,t + T))(A z (x',t)) 
e 2 (n) 



At 



dt'dv P(w)fc(x, t' + t; t>)fc(x', t'; v), (39) 



where the 'pulse function' k(. . . , t; . . .) is defined in Eq.(p7|). 
The prefactor e 2 (n) /At = el is the shot noise spectrum SNj 
(Eq.|l0|). The Fourier integral with respect to the time differ- 
ence r gives the squared Fourier transform of the pulse func- 
tion that can be evaluated analytically. We do not give this 
formula here, but proceed directly to the magnetic field cor- 
relation tensor. After differentiation, one needs the Fourier 
transform 



dt- 



ve 



2e iujz / v \uj/v\ 



r 2 + (z- vt) 2 f' 2 



Ki{r\u/v\), 



(40) 

where r 2 = x 2 + y 2 and K i is a Bessel function. The mag- 
netic noise field has the same orientation (azimuthal) as the 
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static field so that in cylindrical coordinates, the correlation 
tensor has a single nonzero element given by 



S vtp > (x,x'; oj) 



Mo SN X 



An 2 rr 



dvP(v)e i(z - z ' )uj/v 



(41) 



In the low-frequency limit where u>h,(jr <C v, one gets 
(rui/v)Ki(rcu/v) « 1 + 0((ru/v) 2 ). At the guide center 
r = r' = h, z = z', we then recover the magnetic noise 
spectrum ([ll]) with Si(lu) = SNi. 

The scattering cross section (^8j) involves the Fourier 
transform of the correlation function ( |4l| ) with respect to the 
distance s = z — z' . This gives a J-function that permits 
us to perform the integral over the electron velocities. At 
r = r' = h, we get 



dvP 



Ki(h\u/v\) 



g 



(42) 
(43) 



For the scattering process pi — ► ^ + w = — Wfj is the 
negative kinetic energy difference (see Eq.E6h, so that 



g 



q Pi + hq 2 /2 
Mq 



EL 

M' 



(44) 



neglecting the recoil shift compared to the Doppler shift. This 
approximation yields the scattering probability (34). 
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Fig. 1 Principle of a linear magnetic quadrupole guide ("side 
guide"). Figure courtesy of Sierk Potting. 
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Fig. 3 Rate of spin flips in a microtrap above a planar metal sub- 
strate. The black solid lines describe the near field noise gener- 
ated by thermally excited currents (Johnson noise) in the substrate, 
as given by Eq.(^). The dots are an exact numerical calculation 
from [|lj]. Larmor frequencies lul = /xs|B(r)|/?i = 1 MHz x 2tt 
(a) and 100 MHz x 2ir (b) are chosen. The red dash-dotted lines 
describe the noise due to fluctuations in the electric current of a side 
guide. The current noise is assumed at shot noise level (0 dB) or 
20 dB above shot noise. The guide height h is lowered by ramping 
down the wire current / with a constant ratio I/h and at fixed bias 
field B b = 100 G, cf. EqM. 



y=5, 1 =0.32 



Fig. 2 Magnetic near field noise spectrum. The near field spectrum 
is calculated along the lines given in JT^]. The source is a copper 
half-space at T = 300 K, observation distance h — 1 ^im. The trace 
of the magnetic correlation tensor is shown. The top labels give the 
wavelength A. 
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Fig. 4 Spatially averaged coherence function ( |30j ) of a condensate 
expanding in a noisy waveguide. Upper row: with noise giving a 
total scattering rate 7 = 10 and correlation length l c — 1/y/W. 
Lower row: no noise (7 = 0). Left column: noninteracting cloud 
(g — 0), right column: with interactions (g = 10). The dotted lines 
in the upper left panel are the analytical prediction (^l|). The units 
are harmonic oscillator units relative to the initial confinement along 
the guide axis. 
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Fig. 5 Model for shot noise: ballistic electron flow. 
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Fig. 6 Magnetic noise spectrum (^Tj) from the ballistic electron 
model. The spectrum is taken at the side guide center (height h = 
1 fj,m) and normalized to its low-frequency limit. The electron ve- 
locity distribution P(v) is that of Maxwell or Fermi-Dirac with 
characteristic velocities as given in the figure. The distribution is 
centered around the drift velocity vd ■ 



